The Convergence Depth of the Local Pecuhar Velocity 
Field 



00 



Daniel A. Dale 

Infrared Processing and Analysis Center, California Institute of 
Technology, Pasadena, CA 91024, USA 

Riccardo Giovanelli 



W)' Department of Astronomy, Cornell University, Ithaca, NY 14853, USA 

■ Abstract. We have obtained Tully-Fisher (TF) measurements for some 

3000 late-type galaxies in the field and in 76 clusters distributed through- 
^ ■ out the sky between ~10 and 200/i~^ Mpc. The cluster data are applied 

to the construction of an / band TF template, resulting in a relation 
J> ■ with a scatter of ~0.35 magnitudes and a zero-point accurate to 0.02 

magnitudes. Peculiar motions are computed by referral to the template 
relation, and the distribution of line-of-sight cluster peculiar motions is 
presented. 

' The dipole of the refiex motion of the Local Group of galaxies with 

0^ . respect to galaxies with measured peculiar velocity converges to the CMB 

I dipole within less than 6000 km s^^. The progression of this convergence 

• is well illustrated when the refiex motion is referred to a well-distributed 

C^l sample of field galaxies, and it is maintained when the refiex motion is 

Q ■ referred to the reference frame constituted by the distant clusters in our 

pj I sample. The field and cluster samples exhibit bulk motion amplitudes of 

c/? ■ order 200 km s~^ or smaller. 

Finally, we apply our sample of cluster peculiar velocities to a test of 
the putative Hubble bubble recently claimed by Zehavi and coworkers. In 
^ ', contrast to their findings, our data support a relatively quiescent Hubble 

' fiow beyond ~35/i~^ Mpc. 



1. The SFI, SCI and SCII Surveys 

The results discussed in this report derive from peculiar velocity measurements 
obtained for three samples of spiral galaxies, based on the Tully-Fisher (TF) 
technique which combines / band photometry and either 21 cm or Ha long-slit 
spectroscopy. The three samples are: (a) SFI, which includes 1631 field galaxies 
out to cz ~ 6500 km s"^ (Haynes et al. 1998a,b); (b) SCI, which includes 782 
galaxies in the fields of 24 clusters within cz ~ 9000 km s~^ (Giovanelli et al. 
1997a,b); and (c) SCII, which includes 522 galaxies in the fields of 52 clusters 
between cz ~ 5,000 and 20,000 km s"! (Dale et al. 1997, 1998, 1999b,c). For 
the merged cluster sample SCI-I-SCII, there is an overlap of: 49 clusters with 
the Lauer &: Postman (1994) cluster sample; 27 clusters with the SMAC sample 



1 



(Hudson et al. 1999); and 22 clusters with the EFAR sample (Wegner et al. 
1996). 

Detailed criteria for the construction of the samples are given elsewhere: 
for SCI see Giovanelli et al. (1997a); sample characteristics of SCI are given 
in Giovanelli et al. (1994); for SCII see Dale et al. (1999c). Driving concerns 
for the selection of each sample were: dense sampling of the local velocity field 
(SFI), high quality determination of the TF template relation (SCI) and reliable 
recovery of the dipole reflex motion of the Local Group (SCI+SCII). The latter 
sample was constructed with the aid of Monte Carlo realizations of the cluster 
distribution and velocity field, utilizing mock samples from A^-body simulations 
obtained under a variety of cosmological scenarios, with the collaboration of 
Stefano Borgani. 

2. The / Band TuUy-Fisher Template Relation 

Tully-Fisher distances, cztf, and therefore peculiar velocities 

Vpec = cz- CZtf , (1) 

are computed by reference with a fiducial template relation which must be ob- 
servationally derived. The template relation defines the rest reference frame, 
against which peculiar velocities are to be measured. The importance of accu- 
rately calibrating such a tool cannot be overemphasized. It is possible that the 
discrepancies between different claims of bulk motions may be partly related to 
insufficiently well determined template relations. 

For an assumed linear TF relation we need to determine two main param- 
eters: a slope and a magnitude zero point. The slope of the TF relation is best 
determined by a sample that maximizes the dynamic range in absolute mag- 
nitude Mj and disk rotational velocity width i.e. one that preferentially 
includes nearby objects. The proximity of the SCI sample provides the dynamic 
range in galactic properties necessary to accurately determine the TF slope. On 
the other hand, the magnitude zero point of the relation is best obtained from a 
sample of distant objects for which a peculiar motion of given amplitude trans- 
lates into a small magnitude shift. The SCII sample provides an ideal data set 
from which to calibrate the TF zero point. 

2.1. The Accuracy of the Tully-Fisher Zero-Point 

As discussed in Giovanelli et al. (1997b) and Giovanelli et al. (1999), given 
a number A^ of clusters the uncertainty on the TF zero point of the resulting 
template cannot be depressed indefinitely by increasing the average number n 
of galaxies observed per cluster, and taking advantage of the n~^/^ statistical 
reduction of noise on the mean. That is because a "kinematical" or "thermal" 
component of the uncertainty depends on the number A^, the distribution in 
the sky and the peculiar velocity distribution function of the clusters used. In 
SCI, for example, the statistical uncertainty deriving from the total number of 
galaxies observed {n x A") is exceeded by the kinematic uncertainty, which is 

quantified as follows. For a sample of N objects with an rms velocity of (F^)^^^ 
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at a mean redshift of {cz) , the expected accuracy of the zero point is hmited by 
systematic concerns to 



|Am| 



2.17 (y^)'/' 

{cz)VN 



mag. 



(2) 



This quantity is about 0.04 mag for SCI, while it is only 0.01 mag for SCII 

due to the larger mean distance and number of clusters of the latter; the SCII 
sample was selected, in part, to improve upon the accuracy of the kinematical 
component of the TF zero point. Since the total number of galaxies involved in 
the two samples is comparable, the zero point of the SCII template is thus more 
accurate than that of SCI. The total uncertainty of the SCII zero point, including 
the contribution from limitations on its statistical, or internal, accuracy (arising 
from measurement errors, uncertainties in the corrections applied to observed 
parameters, the TF scatter, etc) is of order 0.02 magnitudes. 

The TF template relation is determined internally for a cluster sample. In 
the case of SCI, it was obtained by assuming that the subset of clusters farther 
than 40/i~^ Mpc has a globally null monopole (Giovanelli et al. 1997b). The 
SCII-based template is obtained by assuming that the overall set of clusters 
(wich extends between 50h~^ and 200/i^^ Mpc) has a globally null monopole, 
and adopting the same TF slope as for the SCI sample (Dale et al. 1999c). This 
approach does not, however, affect the value of the dipole or higher moments, and 
thus still allows an effective measure of possible bulk flows. The SCII template 
relation is 



The zero points of the SCI and SCII templates were found to agree to within 
the estimated uncertainty of 0.02 mag (the SCII zero point is fainter by 0.015 
mag). 

2.2. The TuUy-Fisher Scatter and its Intrinsic Component 

Any relation involving observed parameters has a limited accuracy described by 
the amplitude of the relation's scatter. Claims of the scatter in the TF relation 
vary from as low as 0.10 mag (Bernstein et al. 1994) to as high as 0.7 mag 
(Sandage et al. 1994a,b; 1995; Marinoni et al. 1998). The amplitude of the 
scatter does depend on wavelength, and studies in the / band typically yield 
the tightest relations. The efforts with the largest samples yield la dispersion 
values of ~0.3-0.4 mag (Mathewson, Ford & Buckhorn 1994; Willick et al. 1995; 
Giovanelli et al. 1997b; Dale et al. 1999c). Uncertainties in observational mea- 
surements are not the only factors that lead to the overall spread in the data. 
The corrections typically applied to the observed fluxes and disk rotational ve- 
locities are not exactly known, nor are the variety of methods used to account 
for inherent sample biases such as cluster incompleteness. Moreover, there is an 
intrinsic component to the TF dispersion since individual galaxies have diverse 
formation histories. In fact, Eiscnstcin & Loeb (1996) advocate an intrinsic scat- 
ter of 0.3 magnitudes, a number greater than most estimates from observational 
work. In light of this fact, they make the interesting claim that either spirals 
formed quite early or that there must be a type of feedback loop that promotes 
galactic assimilation. The results of Eisenstein &; Loeb are not corroborated by 



Mi - 5 log h 



7.68(log W - 2.5) - 20.91 mag. 



(3) 
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the results of A^-body simulations (e.g. Baugh et al. 1997, Mo et al. 1998 and 
Steinmetz & Navarro 1999), in which reasonably low values for the TF scatter 
are recovered. 

As already established in Giovanelli et al. (1997b), Dale et al. (1999c), and 
Willick (1999), Figure |I| reinforces the notion of low intrinsic scatter. The left 
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Figure 1. The observed scatter in the TF relation for both the SCI 
and sen datasets. Running averages (in units of magnitudes) of the 
observed errors are given as dotted lines. The thin solid line in each 
panel is the sum in quadrature of the two dotted lines and represents 
the total uncertainty resulting from measurement errors and that deriv- 
ing from extinction, geometric and statistical corrections. The circles 
plotted indicate the standard deviations of residuals from the template 
TF relation. The thick solid line in each panel is the total scatter in 
the TF relation and has been computed by a sum in quadrature of the 
observed error distribution and an intrinsic component equivalent to 
that found in Giovanelli et al. (1997b; see Equation ^ below). 

panel refers to the SCI data whereas the righthand panel shows the results from 
sen. The two dotted lines in each panel indicate the velocity width and mag- 
nitude uncertainties Cx and €y, with multiplied by the TF slope b to put it on 

a magnitude scale; the thin solid line labeled = ^ (bex)^ + represents the 
average uncertainty resulting from measurement errors and that deriving from 
extinction, geometric and statistical corrections. The data displayed in Figure 
are generated using equal numbers of galaxies per data point (the difference 
in the redshift distributions of the SCI and SCII samples result in slightly dif- 
ferent observed velocity width distributions). The circles plotted represent the 
average standard deviations of the residuals from the fiducial TF relation. We 
see that the velocity width errors dominate those from the / band fluxes, which 
are approximately independent of velocity width. Furthermore, the logarithmic 
velocity widths become increasingly uncertain for slower rotators (cf. Giovanelli 
et al. 1997b; Willick et al. 1997). We approximate the total observed scatter 
for the SCI and SCII samples with simple linear relations that depend on the 
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velocity width: 



o"tot,sci = -0.33x + 0.32 mag. 
o"tot,scii = -0.40x + 0.38 mag. 



(4) 
(5) 



The total scatter for SCII is in general larger than that for SCI. This is 
unsurprising given that SCII velocity widths primarily stem from optical rotation 
curves instead of 21 cm profiles, a comparatively easier source from which to 
estimate velocity widths, and since the nearer SCI galaxies generally have better 
determined disk inclinations. The gap between the observed scatter and the 
measured errors is attributed to an intrinsic scatter contribution: the thick top 
line is a sum in quadrature of our observed measurement errors, 6^, and an 
intrinsic scatter term (Giovanelli et al. 1997b): 



the same in SCI and SCII. This result has been recently confirmed by Willick 
and his collaborators, as reported at this meeting. 

3. Peculiar Velocity Distributions 

We interpret departures from the canonical template relation as an indication 
of peculiar motion, with larger departures from the template implying larger 
amplitude peculiar velocities. Quantitatively, for an object at a redshift z with 
an average magnitude departure from the template of Am, we write the peculiar 
velocity as 



Tabulations of the SCI and SCII peculiar motions are given in Giovanelli et 
al. (1999b) and Dale et a. (1999c), respectively. Figure ^ shows the distribu- 
tion of SCI and SCII (CMB frame) peculiar velocities in an Aitoff projection 
of Galactic Coordinates. The symbols plotted in the figure reflect both the ra- 
dial directions of the peculiar velocities and the strengths of the measurements 
- in the CMB reference frame, solid/dotted-lined squares represent approach- 
ing/receding clusters and the square size is inversely proportional to the accuracy 
of the measurement. The largest cluster peculiar velocities, e.g. those for A3266 
and A3667, are also the most uncertain, as the clusters are poorly sampled. 

3.1. The RMS One Dimensional Peculiar Velocity Dispersion 

It is useful to estimate the line-of-sight distribution of peculiar velocities. The 
amplitude of that distribution has been known to be a very sensitive discrimi- 
nator of cosmological models (see, for example, Bahcall &: Oh 1996). The SCII 
cluster sample is relatively distant and the cluster membership counts are rela- 
tively anemic when compared to SCI. Consequently, SCII peculiar velocities are 
much less certain and the overall distribution shown in Figure |^ is significantly 
broadened by measurement errors. The peculiar velocities are represented by 
equal area Gaussians centered at the peculiar velocity of each cluster with dis- 
persions equal to the estimated peculiar velocity errors. The thick dashed line 
superimposed on each plot is the sum of the individual Gaussians (its amplitude 



(Tint = -0.28x 0.26 mag. 



(6) 




(7) 
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Figure 2. The all-sky distribution of the SCI and SCII peculiar ve- 
locity samples in Galactic coordinates. The symbol sizes are inversely 
proportional to the peculiar velocity uncertainties. The examples at 
the bottom give the scales. Solid/dotted-lined squares refer to pos- 
itive/negative peculiar velocities in the CMB frame. Asterisks mark 
the apices of the motion of the Local Group with respect to the CMB 
and the Lauer & Postman (1994) cluster inertial frame. 



has been rescaled for plotting purposes). The la dispersion in the observed dis- 
tribution of peculiar velocities is found from a Gaussian fit to the dashed line: 
cid,obs,SCl = 325 km s^^ and crid,obs,SCll = 796 km s^-*^. These values, however, 
are biased high by measurement errors. Recovering an estimate of the true values 
can easily be obtained via Monte Carlo simulations, yielding (Tid,sci = 270 it 54 
km s~^ and o"id,scil = 341 it 93 km s~^ where the error estimates derive from 
the scatter in the dispersions of the simulated samples. These values of aid 
are consistent with a relatively low density Universe (Giovanelli et al. 1998b; 
Bahcall & Oh 1996; Borgani et al. 1997; Watkins 1999; Bahcall, Gramann & 
Cen 1994a; Croft & Efstathiou 1994). 



4. Dipoles of the Field Spiral Sample (SFI) 

On a grid of equatorial coordinates in an Aitoff projection centered at a = 6^* 
and 5 = 0°, Figure |^ displays galaxies from the SFI sample. Filled and unfilled 
symbols refer respectively to positive and negative peculiar velocities, measured 
with respect to the Local Group reference frame. Lines of Galactic latitude 0°, 
+20° and —20° are also shown, outlining the Zone of Avoidance. This sample, 
which extends to 6500 km s~^, clearly displays the dipole moment associated 
with the motion of the Local Group with respect to the CMB, in the form of 
a prevalence of positive peculiar velocities in the Southern Galactic Hemisphere 
and of negative ones in the Northern half (the apex of the Local Group motion 
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Figure 3. A display of equal area Gaussians representing the pecu- 
liar velocity samples for SCI and SCII. Each Gaussian is centered at 
the value of a cluster peculiar velocity, has a dispersion given by the 
uncertainty in the peculiar velocity measurement, and has an ampli- 
tude that is inversely proportional to the uncertainty. The sum of the 
Gaussian profiles is given by the thick dashed line. 

is near a = 11.2^, b = —28°). In other words, a large fraction of the volume 
sampled by the SFI galaxies does not share the Local Group motion resulting 
in the CMB dipole. Figure |^ is a similar display to that in Figure except 
that the peculiar velocities are referred to the CMB reference frame. The dipole 
signature clear in Figure ^ is gone, indicating that most of the SFI sample has a 
small bulk flow with respect to the CMB reference frame. This issue is tackled 
in a more quantitative way in what follows. 

Figure ^ displays the dipoles of the reflex motion of the Local Group with 
respect to field galaxies in shells 2000 km s~^ thick. The dashed line in panel 
^ corresponds to the amplitude of the CMB dipole. The three sets of sym- 
bols identify different ways of computing the peculiar velocities, using different 
subsets of the data or adopting a direct (stars) or inverse Tully-Fisher relation. 
The SFI sample indicates convergence to the CMB dipole, in the motion of the 
Local Group with respect to galaxies in shells, when a radius of a few thou- 
sand km s~^ is reached. Convergence is achieved both in amplitude and apex 
direction. Note that the reflex motion of the LG is estimated with respect to 
shells centered on each redshift, and not with respect to all galaxies within that 
redshift; the convergence rate in the latter case would be slightly slower than 
shown in Figure |^. The dipole of Lauer & Postman (1994) is excluded with a 
high level of confidence. Alternatively, the bulk flow with respect to the CMB 
reference frame of a sphere of 6500 km s~^ is 200ib65 km s~^, directed toward 
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Figure 4. The distribution of SFI peculiar velocities in the Local 
Group reference frame. 
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Figure 5. The distribution of SFI peculiar velocities in the CMB 
reference frame. 
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Figure 6. Amplitudes (a) and apex directions (b) of the peculiar 
velocity dipole for the SFl sample. 'CMB' and 'LP' identify the apices 
of the CMB and Lauer & Postman (1994) dipoles. The dashed line in 
(a) is the amplitude of the CMB dipole. 



{l,b) = (295°, +25°) lb 20°. This is in general agreement with the direction of 
bulk flows reported in other studies (da Costa et al. 1996; Courteau et al. 1993; 
Dekel 1994; Hudson et al. 1999; Wilhck 1999), but it is smaller than other 
determinations, which range between 270 and 700 km s~^. See Giovanelli et 
al. (1998a) for further details. Compare these results with the convergence ex- 
pectations from the PSCz redshift survey, as presented by W. Saunders at this 
meeting. 



5. Dipoles of the Cluster Spiral Samples (SCI and SCII) 

The dipole moments of the reflex motion of the Local Group with respect to (i) 
the subset of SCI clusters farther than 3000 km s^^ and (ii) the SCll sample plus 
the subset of SCI clusters farther than 4500 km s^^ both coincide, within the 
errors, with that of the CMB. This coincidence is matched both in amplitude 
and apex direction. Figure |^ displays the error clouds for the coordinates of the 
dipole derived from the distant SCl+SCll sample, plotted against supergalactic 
Cartesian coordinates. One and two-sigma confidence ellipses are plotted, as well 
as the locations of the CMB and the Lauer &; Postman (1994) reflex dipoles. The 
latter can be excluded by our data with a high level of confidence (at the 3.5a" 
level). The bulk motion for the 64 clusters comprising the distant SCl+SCll 
sample {{cz) ~ lOO/i^^ Mpc) is estimated to be no greater than 200 km s^^; 
corrected for error biasing, the amplitude is consistent with a null bulk flow. In 
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Figure 7. Error clouds of the dipole of the reflex motion of the Lo- 
cal Group with respect to the distant SCI-I-SCII sample. The error 
clouds (with 1-cr and 2-a confidence ellipses) are plotted in supergalac- 
tic Cartesian coordinates. The CMB dipole is a filled symbol, which 
appears enclosed, in all plots, within the 1-a contours; the large cross is 
the dipole of Lauer & Postman (1994), which is excluded by the data. 



other words, the inertial frame defined by the SCI-I-SCII clusters is consistent 
with being at rest with respect to the CMB fiducial rest frame. More locally, the 
bulk flow within a sphere of 6000 km s~^ radius, as derived from SCI, is between 
140 and 320 km s^^, in the CMB reference frame. See Dale et al. (1999a) and 
Giovanelli et al. (1998b) for more details. 



6. No Hubble Bubble in the Local Universe 

Recently, it has been suggested by Zehavi et al. (1998) that the volume within 
cz ~ 7000 km s~^ is subject to an acceleration, in the sense that the local Hub- 
ble constant is 6.6 it 2.2% larger than the fiducial value. The inference is that 
we reside at the center of a local underdensity of amplitude 20%, surrounded by 
an overdense shell. The result is based on the distances to 44 SN of type la. In 
Figure |8|, we test that result by using the combined sample of all 76 SCI-I-SCII 
clusters. The TF distances to each individual cluster have a typical accuracy 
similar to or better than that quoted for the SN measurements (which is 5-8% 
due to internal errors alone), and the cluster sample is well distributed over the 
sky (see Figure ||). Figure § displays 5H/H = Vpec/cztf against hd = cztf/100, 
where Vpec is the peculiar velocity in the CMB reference frame and cztf is the 
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Figure 8. No Hubble Bubble in the local Universe. 5H/H = line 
identifies null deviation from Hubble flow. The dotted line is the effect 
of the Hubble bubble proposed by Zehavi et al. (1998). The lower 
panel shows the error bars associated with the filled data points shown 
above. 



Tully-Fisher distance in Mpc, in the same frame. Unfilled symbols represent 
clusters with poor distance determinations, based on fewer than 5 galaxies with 
Tully-Fisher measurements per cluster. Figure ^ illustrates how, at small dis- 
tances, the deviations from Hubble flow are dominated by the motions of nearby 
groups, comparable in amplitude to those of the Local Group (611 km s~^). At 
small distances, even modest peculiar velocities constitute a sizable fraction of 
cz, thus amplifying and — if the sampling is sparse — even distorting the dis- 
tribution of values of SHo/Hq. The deviation from Hubble flow that they imply 
is of scarce interest. At distances larger than 35/i~^ Mpc, the monopole of the 
cluster peculiar velocity field exhibits no significant change of value. Our sample 
is most sensitive to the possible presence of a "step" (as would be produced by 
the Hubble bubble proposed by Zehavi et al.) between 50/i^^ and 130/i^^ Mpc. 
The amplitude and significance of a step at 70/i^^ Mpc is 

= 0.010 ± 0.022 (8) 

Ho 

Our data yield no evidence for a step deviation in the Hubble flow of amplitude 
larger than 2-3%, up to 130/i~^ Mpc. The Hubble bubble suggestion of Zehavi 
et al. (1998) is thus not corroborated by this larger and more accurate data set 
(see Giovanelli et al. 1999 for further details). 



7. Conclusions 



We have obtained / band Tully-Fisher data for 3000 spiral galaxies in the field 
and in 76 clusters, evenly spread across the sky and distributed between ^-10 and 
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200/t~ Mpc. Rotational velocity widths derive from either 21 cm or long-slit 
optical spectroscopy. The SCI and SCII cluster data are used to construct an 
accurate Tully-Fisher template relation. The relation has an average scatter of 
approximately 0.35 magnitudes, and a zero point with an accuracy of 0.02 mag. 

Peculiar velocities are obtained for the SFI galaxies and each of the SCI 
and SCII clusters, by reference to the TF template relation. The typical dis- 
tance uncertainty is 16% for individual galaxies and 3-9% for clusters, the latter 
primarily depending on the number of Tully-Fisher measurements within each 
cluster. The rms line-of-sight component of the cluster peculiar velocity, de- 
broadened for measurement errors, is 270 ± 54 km for SCI and 341 it 93 km 
s-i for SCII. 

Using the data for the field SFI and cluster SCI samples of spiral galaxies, 
we have shown that the reflex motion of the Local Group converges to the CMB 
dipole amplitude and direction within 6000 km s~^. Results from the relatively 
distant SCII cluster sample confirm that convergence is maintained beyond the 
limits of the SCI and SFI samples. 

Finally, no evidence is found for the local Hubble bubble advocated by 
Zehavi and coworkers, or for any radially averaged deviation from Hubble flow 
with amplitude larger than 2-3%, between 35h~^ and 130h~^ Mpc. Averaged 
over distance in excess of a few tens of Mpc, the Hubble flow appears to be 
remarkably smooth. 
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